arXiv:1503.02165v4 [gr-qc] 20 Jan 2016 


Noname manuscript No. 

(will be inserted by the editor) 


Gravitational Collapse, Shear-free Anisotropic 
Radiating Star 


B.C. Tewari • Kali Char an 


the date of receipt and acceptance should be inserted later 


Abstract We here present a new model of a radiating star for a shear-free 
spherically symmetric anisotropic fluid undergoing radial heat flow collapsing 
under its own gravity. The interior metric fulfilled all the relevant physical 
and thermodynamic conditions and matched with Vaidya exterior metric over 
the boundary. Initially the interior solutions represent a static configuration of 
perfect fluid which then gradually starts evolving into radiating collapse. We 
have observed that the model is well behaved for a set of model parameters 
and there are a number of such sets for which the model is well behaved. The 
apparent luminosity as observed by the distant observer at rest at infinity and 
the effective surface temperature are zero in remote past at the instant when 
the collapse begins and at the stage when collapsing configuration reaches the 
horizon of the black hole. 

Keywords Exact solutions • Anisotropic radiating star • Gravitational 
collapse • Black hole 


1 Introduction 

The study of gravitational collapse attracts researchers towards it by the fact 
that it represents theoretical and observable phenomena in the universe. It 
is an important open issue in relativistic astrophysics whether the final out¬ 
come of the gravitational collapse is formation of a black hole or a naked 
singularity(Joshi-Malafarina [1] and references therein). Still there is no es¬ 
tablished theory available which can determine the final fate of this collapse. 
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There is no iron-clad evidence that black hole candidates are indeed black 
holes. There is no logic that prevents existence of naked singularities and as 
per Cosmic Censorship Conjecture Penrose [2] himself considers this an open 
question. Understanding the characteristics and features of final fate of a col¬ 
lapsing system is not just important from the theoretical perspective; it has 
tremendous observational consequences as well. Virbhadra with his collabo¬ 
rators showed that Black holes and naked singularities could be observation- 
ally differentiated through their gravitational lensing features (Virbhadra-Ellis 
[3,4]; Virbhadra-Keeton [5]; Virbhadra [6] and also references therein). 

In view of classical gravity, to understand the nature of collapse and phys¬ 
ical behaviour of a collapsing system properly, it is necessary to construct a 
realistic model of the collapsing system. It is very difficult task on account of 
the highly non linear nature of the governing field equations. Many efforts have 
been made for reducing the complexity and various methods for simplifications 
are frequently adopted and the pioneering work of Oppenheimer-Snyder [7] 
was the initial step in this direction when collapse of a highly idealised spher¬ 
ically symmetric dust cloud was studied. Since then several attempts have 
been made to construct realistic models of gravitationally collapsing systems 
to understand the properties and nature of collapsing objects. It gained to 
tremendous momentum when Vaidya [8] presented a solution describing the 
exterior gravitational field of a stellar body with outgoing radiation and the 
modified equations in this case proposed by Misner [9] and Lindquist et al. [10] 
for an adiabatic distribution of matter. 

It is an established fact that gravitational collapse is a highly dissipating 
energy process (Herrera-Santos [11]; Herrera et al. [12], Mitra [13],[14]). How¬ 
ever, the dissipation of energy from collapsing fluid distribution is described in 
two limiting cases. The first case describes the free streaming approximation 
and a number of solutions of radiating fluid ball in this case discussed by Tewari 
[15-18]. Further Pant- Tewari [19] presented a quasar model with all degree of 
suitability. While second one is diffusion approximation and in this case the 
dissipation is modeled by heat flow type vector and in this the model proposed 
by Glass [20] has been extensively studied by Santos [21] for the junction con¬ 
ditions of collapsing spherically symmetric shear-free non-adiabatic fluid with 
radial heat flow. On a similar ground a number of studies have been reported 
by de Oliveira et al. [22]; de Oliveira-Santos [23]; Bonnor et al. [24]; Banerjee 
et al. [25]; Govinder-Govender [26]; Maharaj et al. [27]; Herrera et al.[28- 
30]; Naidu-Govender [31]; Sarwe-Tikekar [32]; Ivanov [33]; Pinheiro-Chan [34]; 
Tewari [35, 36]; Tewari-Charan [37, 38] and also references therein for de¬ 
scribing a collapsing fluid radiating energy. A remarkable work for collapsing 
anisotropic radiating star is due to Herrera -Santos [39] explored the properties 
of anisotropic self-gravitating spheres using the perturbation method. Herrera 
[40] investigated that the local pressure anisotropy is one of the responsible 
factors for inhomogeneities in energy density. On a same ground a noticeable 
work is also due to Tikekar and Sliarma [41]; Nguyen-Pedraza [42]; Nguyen- 
Lingam [43]; Sharif-Abbas [44]; Abbas [45]; Tewari-Charan [46], Ivanov ([47], 
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[48]) and also references therein have been reported with the impact of inho¬ 
mogeneity, anisotropy and various dissipative processes on the evolution. 

The main objective of this work is to shed some light on the final fate 
of gravitational collapse by taking Tewari and Charan [46] solution as seed 
solution. Final fate of our model is the formation of a black hole. By taking a 
set of model parameters we have developed a Supernovae model. The interior 
matter fluid is spherically symmetric, shear-free, anisotropic radiating away 
its energy in the form of radial heat and contracting in size during the process 
of collapse. The interior metric matched with the Vaidya exterior metric [8] 
over the boundary. The paper is organised as: In Section 2 the interior space- 
time and related equations for collapsing system are given. In Section 3 the 
exterior space-time and junction conditions are presented for matching the 
two space-times. Section 4 describes a detailed study of the model in which 
we have obtained expressions of various physical parameters and temperature 
profile of the collapsing body. Finally in Section 5 some concluding remarks 
have been made. 


2 The interior space-time and related equations for collapsing 
system 

The interior metric of a shear-free spherically symmetric fluid is given by 

ds 2 _ = — A 2 (r, t)dt 2 + B 2 {r, t){dr 2 + r 2 (d0 2 + sin 2 8d(j) 2 )} (1) 


The energy-momentum tensor for the matter distribution with anisotropy in 
pressure is 


= (e + p^w^Wv + p t g^u + (p r - Pt)x^x„ + q^w v + (2) 


where e is the energy density of the fluid, p r the radial pressure, pt the tan¬ 
gential pressure, is the four- velocity, q M the radial heat flow vector and x M 
is a unit space like four vector along the radial direction. 

Assuming comoving coordinates, we have uT = . The heat flow vector 

is orthogonal to the velocity vector so that = 0 and takes the form 

q 1 * = q5i -The line element (1) corresponds to shear- free spherically symmetric 
fluid (Glass [49]), as the shear tensor vanishes identically. 

In order to solve the Non-trivial Einsteins field equations which are gen¬ 
erated by (1) and (2), we choose a particular form of the metric coefficients 
given in (1) into functions of r and t coordinates as A(r,t) = Ao(r)g(t) and 
B(r,t) = Bo(r)f(t). The coupling constant in geometrized units is taken as 
k = 8n (i.e.G = c = 1) and in view of (1) and (2) with the help of above metric 
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coefficients we get the following expressions for field equations 
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Here the quantities with the suffix 0 corresponds to the static star model with 
metric components Ao(r ), B$(r) and the primes and dots stand for differen¬ 
tiation with respect to r and t respectively. 

To find a new parametric class of exact solutions of pressure anisotropy 
equation which is created by the equations (4) and (5), Tewari-Charan [46] 
assumed that the anisotropy evolves n{p t — p r ) = Z\(r, t) = and got 

the following time independent differential equation 


K , a? 4, B' 0 2B'J 2A' 0 B' 
A 0 B 0 rA 0 rB 0 B% A 0 B 0 


( 10 ) 


Making an adhoc relationship between the variables in (10), Tewari- Charan 
[46] obtained the following solution 


Ho = C 4 (l + C 3 r 2 )(l + Cir 2 )* 

( 11 ) 

£o = C 2 (l + C' 1 r 2 )m 

( 12 ) 

X( r \ _ (2«-2) 4C 3 Cir 2 

°y > (1+ 1) (l-t-C3r 2 )(l+Cir 2 ) 

(13) 


where n, l, C\, C 2 , C 3 and C 4 are constants and 

n= ^{(1 +3) ±(^ 2 + 10 ^ + 17)5} 
where n is real if l > —5 + 2p2 or l < —5 — 2p2. 


( 14 ) 
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3 The exterior space-time and junction conditions 


The exterior space-time of a collapsing radiating star is described by Vaidya 
exterior metric [8] 


ds\ = - (l - 2M ^ ) dv 2 - 2 dRdv + R 2 {d6 2 + sin 2 ddcf 2 ) (15) 

where v is the retarded time and M(v) is the exterior Vaidya mass. 

The junction conditions for radiating star matching two line elements (1) and 
(15) at the boundary continuously across a spherically symmetric time-like 
hyper surface £ are very well known and given by Santos [21] 

(rB) s = R s (v) = TZ(t) (16) 

{Pr)s = (qB)s (17) 

m s ( r , t ) = M(v) = ( - r 2 B' - r ^f\ (18) 

l E 


where is the mass function calculated in the interior at r = rs (Cahill et 
al. [50]; Misner-Sharp [51]). 

The surface luminosity and the boundary redshift zs observed on £ are ob¬ 
tained by de Oliveira et al. [22] 


Ze 


Le = f {r 2 B 3 q} s 




The total luminosity for an observer at rest at infinity is 

dM _ L e 
~ (1 + ze) 2 


(19) 

( 20 ) 


( 21 ) 


In the absence of non-adiabatic dissipative forces the equation (17), (p r )s = 
( qB)s, reduces to the condition [(p r )o]i; = 0 and yields at r = rs = Re 


where 


2 / / 2 2 gf _ 2 ggf 

f f 2 gf f 2 


( 22 ) 


(23) 


Motivated by de Oliveira et al. [22] and Bonnor et al. [24], Tewari [36] assume 
g{t) = f(t) and got the solution of (22). The solution so obtained is identical 
to the solution presented by de Oliveira et al. [22]; Bonnor et al. [24] with 
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g(t ) = 1. Keeping in mind the collapsing configuration range finally Tewari 
[36] obtained 

/ = -2av?(l - VI) (24) 

t = i/n( 1 - Vf) (25) 

We observed that the function /(f) decreases monotonically from the value 
/(f) = 1 at t = —oo to /(f) = 0 at f = 0. It interpret that the collapse of 
the configuration begins in the remote past and gradually starts evolve into 
radiating collapse. 


4 Detailed study of general model of collapsing radiating star 

Using (3)-(9), (11),(12) and (24)and (25), the general expressions for the en¬ 
ergy density, both the pressures and the heat flux for the collapsing body 
are 
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where eo, (p r )o, (pt)o are energy density, radial and tangential pressures 
respectively in static position of the star and they are given by 
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By utilizing (11), (12) and (23), we have 


« = 2 C4rr(1 c^+i1 ) ~ tT1 -[ nC, i( 1 + ^1) + (l + 1)03(1 + C ir |)] (33) 
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4.1 Detailed study of a specific model 

For different values of n, one can obtained a number of solutions using equa¬ 
tions (11) and (12). For n = 0, we get the homogeneous density and anisotropic 
pressure, for n = 1, we get the homogeneous density and isotropic pressure 
and for n = — 1, we get the isotropic pressure. In this way a number of differ¬ 
ent radiating star models can be develop from the above obtained solution. In 
order to maintain the inhomogeniety and anisotropy a Horizon-free case for 
n = —1 — y/2 has been studied by Tewari-Charan [48]. In the present study we 
assume n = —2, to maintain the inhomogeniety and anisotropy for a collapsing 
radiating star, so utilizing (11)-(14), (30)-(32) we get 


(Pr) 0 — 

(Pt) 0 = 


A 0 = C 4 (l + C 3 r 2 )(l+C 1 r 2 ) 2 T 
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49C 2 (1+Cir 2 )¥ + 4Cir j + 
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fCir 2 )(7+] 
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(34) 

(35) 

(36) 

(37) 

(38) 

(39) 


The junction condition {(p r )o}i; = 0 gives 


-C' 1 (7 + 4(7^1) 

7(1 + C ir %)(7 + 5Cir|) + Cir|(7 + 4Cir|) 


(40) 


Here from (38) and (39), we are seeing that at the centre radial and tangential 
pressures are equal and anisotropy vanishes there. 

A physically reasonable solution should satisfy certain conditions and which 

are: 

The central values of both the pressures, density and gravitational potential 
component should be non-zero positive definite. This condition gets satisfied 
if C\ > 0, C 2 > 0, C 4 >0 and C 3 > {—C\)/7 , subjecting to the condition 
^ < 1 which yields C 3 < 

Secondly the solution should have monotonically decreasing expressions for 
both the pressures and density with the increase of radial coordinate r. In view 
of this we have observed that these conditions hold good if C\ > 0 and C 3 < 0 
satisfying all other requirement of the model. Further, it is mentioned here 
that the boundary of the collapsing radiating star is established only when 
Ci > 0, C 2 > 0, C 4 > 0 and < c 3 < 0 . 
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Now utilizing (26)-(29), (34) and (35) the expressions for e,p r ,p t ,q reduce 
to the following 

12a 2 (l-7/) 2 


ne= p 
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The fluid collapse rate 0 = w? with the help of (1), (24), (34) and (35) is 
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where by using (33) 

a = + C 3(7 + 9Cr4)] (46) 

7C 2 (1 + Cir^)? 

To develop a Supernovae model and pictorial representation of various physical 
parameters keeping in mind the reasonable conditions we have chosen the 
model parameters as C\ = 1, C 2 = 0.1, M 0 = 5 Mq, C 4 = 0.5 



Fig. 1 Behavior of density vs radial and temporal co-ordinates 






















Gravitational Collapse, Shear-free Anisotropic Radiating Star 


9 



Fig. 2 Behavior of radial pressure vs radial and temporal co-ordinates 




Fig. 4 Behavior of heat flux vs radial and temporal co-ordinates 
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Fig. 5 Behavior of collapse rate vs radial and temporal co-ordinates 


We can see [fig. 1-4] the physical parameters are finite, positive, monotonically 
decreasing at any instant with respect to radial coordinate for 0 < r < rs- The 
plots of various physical parameters show that the nature of these parameters 
is well behaved with respect to radial and temporal coordinates for a set of 
constraints and there are a number of such constraints for which the model is 
well behaved. Fig. 5 shows that initially collapse rate is zero and it increases 
with respect to radial and temporal coordinates till the formation of black 
hole. 


4.2 Mass energy, physical radius and time of black hole formation 


The total mass energy entrapped inside the surface £ is given by (18), which 
becomes by using (24), (34), (35) and (46) 


M(v) 


8 + C 3 (7 + 9Cir|)] 2 (l - y/f) 2 
49(1 + C'ir|.) J f i (l + Czr 2 s ) 2 


+ m 0 f 
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(47) 


where 


_ 2C 1 C 2 4(7 + 6C 1 r|) 
19(1 • r ,/•;() 7 


(48) 


Using (16) and (35), we get the physical radius of the collapsing radiating star 
as 


Rs(y) =r s C 2 (l+C 1 rl)~ l ' 7 f 


( 49 ) 
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Fig. 6 Behavior of mass vs temporal co-ordinate 
radius 



Fig. 7 Behavior of radius vs temporal co-ordinate 


Utilizing (19)-(21), (24), (34), (35) and (46) the surface luminosity, the 
boundary redshift on £ and the luminosity for distant observer at rest at 
infinity are 
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The expression (52) shows that vanishes in the beginning when f(t) —► 1 
and at the stage when [zs —> oo]. The black hole formation time when the 
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collapse reaches the horizon of the black hole occurs when the surface redshift 
goes to infinity, for this the term in the parentheses of (51) goes to zero and 
we obtain 


n = _ 4C 1 r|[2C 1 + C 3 (7 + 9Cir|)] _ 

VJbh 4Cir|[2Ci+C 3 (7 + 9Cir|)] + (7 + 5C\r|)(l + C 3 r|) 1 ’ 

and 


tBH 


_ (7 + 5C , ir|)(l + C 3 r%) _ 

a 4Cir|[2Ci + C 3 (7 + 9C ir %)\ + (7 + 5C'ir|)(l + C 3 r%) 


(54) 


Figs. (6 and 7) show the behavior of mass energy and physical radius of col¬ 
lapsing radiating star. Initially when /(f) —> 1 both the mass and physical 
radius of collapsing body are 5 Mq and 20.4 km and at the time of black hole 
formation they remain 2.4 Mq and 4.86 km respectively. In view of (53) the 
time of black hole formation is observed as 0.2399 S. 


4.3 Temperature profile for collapsing radiating star 


To obtain the temperature inside and on the surface, we utilize temperature 
gradient law governing the heat transport within the collapsing matter (Israel 
et al. [52]; Maartens [53]; and Martinez [54]) given as 

t(j^ v + w»w v )w a q v . a + q* = - W + w^w v )[T v + Tw v \ (55) 


where K is the thermal conductivity and r is the relaxation time. 

To get a simple estimate of the temperature evolution, we set relaxation 
time as zero in (55) and get 
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If we assume thermal conductivity K = > 0, where 7 and i? are positive 

constants, then on integration (56) we get 


o+ i_T 0 (t) 4(1? + 1) a 1 — y/J 

Ag +1 jn A 0 /3/a 1 j 

where To(t) is an arbitrary function of t. 

The effective surface temperature observed by external observer can be 
calculated from the following expression given by Schwarzschild [55] as 

Tb = { Tr6{rB 0 f ) 2 I s L °° 

= A_Cfr|_ (1 ~ VJ) 1 

49 7r( JC , |( 1 + C , 1 r| : )¥( 1 + C , 3 r|)2 f 5 / 2 (1 + z s ) 2 


( 58 ) 
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where for photons the constant S is given by 5 = (7r 2 k 4 )/(15 h 3 ), here k and h 
denoting respectively Boltzmann and Plank constants. 

Choosing ft = 3 which represents radiation interaction with matter fluid 
through the diffusive approximation (Misner and Sharp [56]). The arbitrary 
function T 0 (t) is determined by using (57) and (58) as 


To (t) = 


C' 4 (l + C' 3 r 2 )(l + C' 1 r 2 ) 2 / 7 1 (1-vT)' 



(59) 


/ 5/2 J^(l + ^) 2 


Temperature distribution throughout the interior of the collapsing body is 
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(60) 





Fig. 8 Behavior of surface temperature vs temporal co-ordinate 


It follows that the surface temperature of the collapsing star tends to zero at 
the beginning of the collapse [/ —► 1] and at the stage of formation of black 
hole [zs —> oo], it is also evident from the fig. (8). 


5 Discussion and Conclusions 

We have given a new class of exact solutions of spherically symmetric shear-free 
anisotropic fluid distributions radiating away its energy in the form of radial 
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heat flow. We have obtained a class of exact solutions by assigning different 
values to the parameter n. A simple radiating star model for n = —2 studied in 
detail. The model is physically and thermodynamically sound as it corresponds 
to well-behaved nature for the fluid density, both the radial and tangential 
pressures and the radiation heat flux throughout the fluid sphere. Initially 
the interior solutions represent a static configuration of dissipative fluid which 
then gradually starts evolving into radiating collapse. The apparent luminosity 
as observed by the distant observer at rest at infinity is zero in remote past at 
the instance when the collapse begins and at the stage of black hole formation. 

We observed that the function /(f) decreases monotonically from the value 
/(f) = 1 at f = —oo to /(f) = 0 at t = 0. The plots of various physical pa¬ 
rameters [figs. (1-4)] show that the nature of these parameters is well behaved 
with respect to radial and temporal coordinates for a set of constraints and 
there are a number of such sets for which the model is well behaved. Fig. (5) 
shows that initially collapse rate is zero and it increases with respect to ra¬ 
dial and temporal coordinates till the formation of black hole. Figs.(6 and 7) 
show the behavior of mass energy and physical radius of collapsing radiating 
star. Initially when /(f) —► 1 both the mass and physical radius of a core of 
pre Supernovae were 5 Mq and 20.4 km and they remain 2.4 Mq and 4.86 
km respectively at the time of black hole formation. The time of black hole 
formation is observed as 0.2399 S. The surface temperature of the collapsing 
radiating star tends to zero at the beginning of the collapse [/ —> 1] and the 
stage of formation of black hole [zs —> oo]. 
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